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We propose a d-density wave induced by the spin-orbit coupling in the CuO plane. The 
spectral function of high-temperature superconductors in the under doped and lightly doped 
regions is calculated in order to explain the Fermi arc spectra observed recently by angle- 
resolved photoemission spectroscopy. We take into account the tilting of CuO octahedra as 
well as the on-site Coulomb repulsive interaction; the tilted octahedra induce the staggered 
transfer integral between p x , y orbitals and Cu t2 g orbitals, and bring about nontrivial effects of 
spin-orbit coupling for the d electrons in the CuO plane. The spectral weight shows a peak at 
around (tt/2,tt/2) for light doping and extends around this point forming an arc as the carrier 
density increases, where the spectra for light doping grow continuously to be the spectra in 
the optimally doped region. This behavior significantly agrees with that of the angle-resolved 
photoemission spectroscopy spectra. Furthermore, the spin-orbit term and staggered transfer 
effectively induce a flux state, a pseudo-gap with time-reversal symmetry breaking. We have a 
nodal metallic state in the light-doping case since the pseudogap has a d x 2_ y 2 symmetry. 
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In recent years, oxide high-T c superconductors have 
been investigated intensively; anomalous metallic behav- 
iors as well as high critical temperatures T c have been 
focused in the study of high-T c cuprates. Clarifying the 
origin of an anomalous metal with a pseudogap is also 
a challenging problem attracting many physicists. Re- 
cently, a peak crossing the Fermi level in the node direc- 
tion of the d-wave gap has been observed in a lightly 
doped La2~ a; Sr :E Cu04 (LSCO) by angle-resolved pho- 
toemission spectroscopy (ARPES). The spectral weight 
shows a peak at around (tt/2,tt/2) for the light-doping 
case and it extends along the Fermi surface with in- 
creasing carrier density. [1,2] The existence of incommen- 
surate correlations has also been reported by neutron- 
scattering measurements, suggesting vertical stripes in 
an underdoped region and diagonal stripes in a lightly 
doped region. Modulation vectors are Q s = (tt ± 2tt5, tt) , 
Q c = (±47rS,0) (or Q s = (tt,tt± 2tt<5), Q c = (0,±4tt<5)) 
for the hole-doping rate x > 0.05, where 5 indicates the 
approximately linear dependence 5 = a;. [3-5] 

In the lightly doped region for x < 0.05, the suggested 
modulation vector is Q s — (tt ± 27r<5, 7r ± 27r<5) and the 
deviation from the linear dependence exhibits S < x.[4] 
From the experiments of resistivity in the lightly-doped 
region of LSCO, the system holds a metallic behavior 
below Tn , which may be due to the formation of metallic 
charge stripes. [6] 

The purpose of this paper is to investigate the spec- 
tral function of doped high-temperature superconductor 
in order to elucidate the recently observed ARPES Fermi 
arc spectra, taking into account the stripe orders sug- 
gested by neutron-scattering measurements. We empha- 
size the importance of distortions of CuO octahedra. 

The tilting of CuO octahedra induces transfer integrals 
between p XtV orbitals and Cu t2 g orbitals, and then the 



spin-orbit coupling between t2 g orbitals is not decoupled 
from e g networks. [8,9] For the spin-orbit coupling 



H, 



(1) 



matrix elements exist between tig orbitals 



(d xz (r) T \H so \d yz (r) T) = 




(2) 


(d yz (r) T \H so \d xz (v) T) = 


me, 


(3) 


(d xz (r) I \H so \d yz {r) |) = 


m% 


(4) 


{d yz (r) I \H so \d xz (r) [) = 


-w% 


(5) 



We have also matrix elements between ti g orbitals and 
e g orbitals: 

(d xz (r) i \H so \d x2 _ y2 (r) T) = £/2, (6) 

(d yz (r) i \H so \d x 2_ y2 (v) T) = -i£/2, (7) 

and those for reversed spins obtained by multiplying by 
— 1. The matrix elements induced by the tilting are 



{p x (x - a/2,y)a\H pd \d xz (r)a) = -t a 



(8) 

(p y (x, y - a/2)a\H pd \d yz {r)a) = -t yz e^ T , (9) 

where r = [x, y), Q = (tt, tt), a is the lattice constant and 
H p d denotes the hybridization term. The factor e lC *' r in- 
duced by the tilting of octahedra in a staggered manner, 
which leads to the doubled unit cell. t xz and t yz are as- 
sumed to be given by sm6 for the tilt angle 0. In the 
low-temperature tetragonal (LTT) phase, the spin-orbit 
coupling may be significantly smaller than that in the 
low-temperature orthorhombic (LTO) phase, since the 
integrals between oxygen p and Cu t2 g orbitals remain 
zero along the tilt axis, where oxygen atoms never move: 
tyz — if the tilt axis is in the y direction. 
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The dispersion in the presence of distortions for the 
five-band p-d model is shown in Fig.l. The parameters 



are as follows: 



t 



pp 

ipd 



0.2, £ 



0.1. 



>. = -2, 
and t T , 



= 0, 



= t 



yz 



= 0.3 in units of 



leV" where t p d is the transfer between d and p 
orbitals and t pp is that between neighboring oxygen p 
orbitals. For the parameters shown above, the splitting 
at X = (tt,0) is of the order of lOmeVwlOOK. Each 
curve is a twofold degeneracy, i.e. Kramers degeneracy 
as it should be. This structure near the Fermi energy is 
well understood, using the single-band model with the 
reduced Brillouin zone, as 



(10) 



where 



£fc = — 2t(cos(fc x ) + cos(ky)) — At 1 cos ( fox ) cos (fey) , (11) 



and Q denotes the wave vector Q 
complex parameter satisfying 



(ir, 7r). A k(T is a 



= A 



k<r- 



A k(T is taken as[10,ll] 

A k(T = icrsin</> • (-2i)Y(k) 



(12) 



(13) 



with F(k) = cos(fc x ) — cos(fc y ) in order to reproduce the 
dispersion for the five-band CuO model near the Fermi 
energy, where <fr = \9; 6 is the tilt angle between the 
Cu-0 plane and the Cu-0 bond, and A is a constant 
estimated as A w 0.2. [11] The dispersion relation of the 
nonintcracting part H n in cq. (10) is 



^k 



6 



k+Q 



(a-a +Q ) 2 + 4|A k(T |2]. (14) 



At half-filling, the low-energy excitations are described 
by the Dirac fcrmion since there is a Fermi point at 
(±7r/2, ±7r/2). A light hole doping results in a small 
Fermi surface around this point with the excitation gap 
near (ir, 0). This peculiar feature induces a pseudogap 
in the density of states and we obtain a nodal metallic 
state. The pseudogap induced by the spin-orbit coupling 
has a d x 2_ y 2 symmetry as apparent from that of A k(T . 
The density of states for small 4> shown in Fig. 2 clearly 
indicates a pseudogap structure for a small excitation 
energy. 

The eigenstate of H resembles the d-density wave pro- 
posed for an anomalous metallic state in the underdoped 
high-T c cuprates[12], and possesses the following order 
parameter 

(15) 



iA DDW Y(k) = {c\ +Qa c ka ) 



for a real constant Addw- Although the d-density wave 
state appears as a solution of the mean-field equations, 
this state is hardly stabilized in variational Monte Carlo 
calculations. Thus, we are motivated to consider the lat- 
tice distortion, which is going to stabilize the d-density 
wave cooperating with the spin-orbit coupling. 

The inhomogeneous ground state under the lattice 
distortion in the underdoped region has been investi- 
gated intensively using the two-dimensional (2D) Hub- 
bard model. [13-16] Within the 2D Hubbard model, the 




Fig. 1. Dispersion relation for the three-band d-p model. The pa- 
rameters are as follows: 2 2 = —2, e p = 2, td asz = e d % = 



-1, tpp = 0.2, J = 0.1, and t x 



: 0.3 in units of t p d, where 



tpd is the transfer between d and p orbitals and t pp is that be- 
tween neighboring oxygen p orbitals. The splitting at X = (n,0) 
is of the order of lOmeV. Each curve is a twofold degeneracy, i.e., 
Kramers degeneracy. 
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4> = 0.037t 
4> = 0.017C 
j ) = 



Fig. 2. Density of states for = 0, O.Olvr and 0.03vr. We set t' = 0. 



linear dependence of S on x has been explained by vari- 
ational Monte Carlo (VMC) methods[16], and further- 
more, the saturation of incommensurability S for x > 
0.125 is also consistent with those obtained by the VMC 
methods. There is also a tendency towards the formation 
of stripes under the lattice distortions, with vertical or 
horizontal hole-rich arrays coexisting with incommensu- 
rate magnetism and superconductivity (SC).[7] 

We are going to evaluate spectral functions in the 
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lightly to optimally doped regions. The tilt angle for 
LSCO is estimated as 14°-18° by EXAFS.[17] Thus, 
9 is approximately 0.17T. If we use A ~ 0.2, we have 
4> = AO « 0.02-7T. Here, we use a slightly larger value, 
4> = 0.057T, in actual calculations to obtain sufficient pre- 
cision because of the numerical difficulty for the small 
splitting. Clearly, we can expect that this does not change 
the peculiar feature of the spectra. We determine the 
variational parameters g and Aq so as to minimize the 
ground state energy for U = 4 and the carrier density 
in the range of < x < 0.2 by the VMC methods. The 
Hamiltonian is 



H = H o + U riifTiii, 



(16) 



where Hq is given in eq. (1). The wave function is written 
in a Gutzwiller form: 



ip = P G lpMF- 
Pg is the Gutzwiller operator 

Pg = !](! - (i - 9)^nn), 



(17) 



(18) 



and the mean field wave function ipMF is obtained as an 
eigenfunction of the Hartree-Fock Hamiltonian 



H MF = H + ^2[5m -sign(a)(~l) 



-u 



'mjlcLc 



(19) 



Srii and m, are expressed by the modulation vectors 
Q s and Q c for the spin and charge part, respectively. [7] 
Equivalently, we use the form given by [7, 18] 



dn,i = -a ^ cosh((arj - xf r )/£ c ), 



and 



in. 



AgJJtanhCCari-afVC.)- 



(20) 



(21) 



The diagonally striped state has hole arrays in the diag- 
onal direction, while the vertically striped state has hole 
arrays in the direction parallel to the x or y direction. Al- 
ternatively, it is also possible to determine the order pa- 
rameter AtQ s a = Y,k( c k+eQ s a C kv) = 0, 1, • • ■ ,M - 1) 
for 6 = 1/M consistently. [19] We obtain the diagonally 
striped state for x < 0.05 and the vertically striped 
state for x > 0.05 as the ground state by the VMC 
methods. [16] Among the diagonal stripes, the bond- 
centered striped state is the most stable when we vary 
the wave function from the site-centered stripe to the 
bond-centered stripe by the VMC methods. We mention 
here that the spin-orbit coupling stabilizes the diagonally 
striped state in the light doping region. [20] 

In the evaluations of spectral functions, we consider 
the effect of Pg within the mean field theory since we 
must consider the excited states as well as the ground 
state. We then evaluate the spectral weight from the 
eigenvalues E am and eigenfunctions (u am )j of the Hamil- 
tonian (H% IF )ij in the real-space representation, where 
i and j are site indices. Green's functions are 



g a (i,j,iuj) = ^2 ' 



(22) 




Fig. 3. Contour map of density of states for the diagonally striped 
state at doping rate x = 0.03 with <f>/n = 0.05. Aq = 0.08, a = 
and t' = -0.2. 
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Fig. 4. Density of states at doping rate x = 0.0612 with <j>/n = 
0.05. Vertical stripes with a 16-lattice periodicity are assumed. 
Aq = 0.10, a = and f = -0.2. 




Fig. 5. Density of states at doping rate x = 0.125 with (j>/ir = 
0.05. Vertical stripes with an 8-lattice periodicity are assumed. 
Aq = 0.16, o = and t' = -0.2. 
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Fig. 6. Density of states for the vertically striped state at doping 
rate x = 0.197 with <j>/n = 0.05. Vertical stripes with an 8-lattice 
periodicity are assumed. Aq = 0.08, a = and t' = —0.2. 
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Fig. 7. Density of states for the normal state at doping rate x = 
0.197 with cfe/n = 0.05, Aq = and i' = -0.2. 

5ff (k,zc) = -L^ e - lk -( R *~ R ^ CT (z,^), (23) 

y 

where ./V a is the number of atoms. The spectral function 
is calculated using the formula 

N„(k,e) =--Img a (k,e + iS). (24) 

7T 

We show the results for x = 0.03, 0.061, 0.125 and 
0.197 for the parameters obtained by the VMC meth- 
ods. The contour map of spectra for the light-doping 
x = 0.03 is shown in Fig. 3 where the calculations 
were performed on a 60 x 60 lattice for the half-filled 
diagonal stripes. A peak near (7r/2,7r/2) appears due 
to the spin-orbit and distortion effects as presented in 
Fig. 3 for — 0.05, while it should be noted that the 
spectra of diagonal stripes without distortion exhibit a 
one-dimensional structure in the diagonal direction. It 
has been pointed out that the diagonally striped state 
can alone explain the opening of the gap if we consider 
the bond-centered stripe. [21] The spectral functions for 
x = 0.061 and x — 0.125 with <fi/ir = 0.05 are shown 
in Figs. 4 and 5, respectively, where we have vertically 
striped states that have 16-lattice and 8-lattice periodic- 
ities, respectively, in accordance with neutron scattering 
measurements. [4] In Fig. 6 the spectral map at x = 0.197 
in the overdoped region is shown, where the stripes have 
an 8-lattice periodicity. We show the spectra without 
stripes at x = 0.197 in Fig. 7 for comparison. We observe 
the absence of spectral weight near (it, 0), which is a char- 
acteristic structure originating from the pseudogap. The 
vertically striped state has one-dimensional-like spectra 
near (±7r,0) with the Fermi wave number kp correspond- 
ing to the one-dimensional quarter-filled band, [13] while 
the Afc term contributes to a peak structure at around 
(±7r/2, ±7r/2). As a result, we obtain the arclike spectra 
for vertical stripes. Thus, as the doping rate x increases, 
the spectra near (ir/2,n/2) in the light-doping case ex- 
tends towards the 2D-like Fermi surface in the optimally 
doped region, which occurs as a crossover. 

In this paper we have examined novel phenomena 
stemming from the spin-orbit coupling induced by the 
tilting of CuO octahedra. We have shown that the char- 
acteristics of the spectral function of doped high-T c 
cuprates can be consistently explained using the 2D elec- 
tronic model with lattice distortions, which is in con- 



5 trast to a theory considering p-orbitals in apical oxygen 
4 atoms. [22] We have taken into account the incommen- 
3 surate structure observed by neutron-scattering exper- 
2 iments. Furthermore, we can expect the following fas- 
-| cinating physics: pseudogap, nodal metal, time-reversal 
symmetry breaking, diagonal stripes, and string-density 
wave as a generalization of the <i-density wave. In the 
half-filled case, we obtain a Fermi point, and thus, a peak 
exists near (tt/2,tt/2) in the density of states that ex- 
tends to form the Fermi arc spectra as the doping rate 
increases. [2] In the low-doping case, we obtain a nodal 
metallic state in the diagonally striped state, since the 
pseudogap has a d x 2_ y 2 symmetry, which is consistent 
with the experiments of resistivity. [6] The arc is expected 
to expand with increasing temperature to form the full 
Fermi surface above the splitting gap ~ 100K.[1] Lastly, 
we comment on a possibility of superconductivity along 
stripes. The coexistence of superconductivity and stripes 
has been pointed out for vertical stripes by the VMC 
methods. [7] It is difficult to have a stable coexistent state 
of SC and stripes for diagonal stripes in the VMC meth- 
ods because the SC pairs must have a d xy symmetry 
along stripes. Thus, superconductivity is suppressed for 
light doping. 

We are grateful to H. Eisaki for valuable discussions. 
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